Empirical Solar Wind Forecasting from the Chromosphere 
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Recently, we correlated the inferred structure of the solar chromospheric plasma 
topography with solar wind velocity and composition data measured at 1 AU. We 
q . now offer a physical justification of these relationships and present initial results of a 

£j , empirical prediction model based on them. While still limited by the fundamentally 

complex physics behind the origins of the solar wind and how its structure develops 
in the magnetic photosphere and expands into the heliosphere, our model provides a 
near continuous range of solar wind speeds and composition quantities that are simply 
estimated from the inferred structure of the chromosphere. We suggest that the derived 
quantities may provide input to other, more sophisticated, prediction tools or models 
such as those to study Coronal Mass Ejections (CME) propagation and Solar Energetic 
Particle (SEP) generation. 

Subject headings: Sun: chromosphere - Sun: solar wind - Sun: solar-terrestrial rela- 
tions 



1. Introduction 



Several authors have attempted to derive simple relations co nnecting large scale vari ations in 
observed solar quantities and plasma quantities measured in situ: iWang & Sheeleyl (U990) showed 
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that solar wind speeds are anti-correla ted with the super-rad ial expansion factors of magnetic flux 
tubes rooted in the solar photosphere; lArge & Piz zo (2000) deve l oped a near-real prediction time 
tool based on lWang & Sheeleyr s derivation. ISchrijver & DeRosal (|2003|) developed si milar predic 



tive t ools using observations from the Solar and Heliospheric Observatory (SOHO; Fleck et al 



1995) . While it is widely accepted that the steady fast and slow solar wind pattern is modulated 
by the rotation of large open and closed magnetic structures on the solar surface, it is clear that the 
solar wind structure is permeated by much small scale structure. 

As yet, it is unknown how these small spatial structures in the solar wind originate; however, 
one reasonable explanation may lie in the fact that these large magnetic structures on the Sun 
are partitioned into even smaller magnetic structures which exist on spatial scales that are below 
current instrumental resolution. If the smaller structures that are embedded in the slow-fast wind 
are the se eds of SEP eve nts driven by large disturbances of the wind (i.e., Coronal Mass Ejections 
or CME s ; iReamesI 1 1 999b then we need to explain, or at least quantify, the continuum of variation 
observed in solar wind velocities and compositions before we can reliably predict the onset of such 
destructive events. 

We present the first results of an empirical model for prediction of solar wind conditions that 
is based on detailed measure ments of chrom ospheric st ructure. This simple model w as developed 
to test the measurements of [Mcintosh et al.1 (|2004|) and iMcIntosh & Leamonl (|2005|) and examine 
why the chromosp heric plasma observed by the Transition Region and Coronal Explorer (TRACE; 
Handy et al.lll999|) appears to couple so well t o that of the helio sphere measured in situ at 1 AU 
by the Advanced Composition Explorer (ACE; IStone et al.lll998|) . We will show that the empirical 
forward propagation of chromospheric structure, through a series of derived power-laws, allows us 
to recover a broad spectrum of variation in solar wind speed and composition which qualitatively 
reproduces that observed. 



2. The Chromospheric Footprint of the Solar Wind 



Mcintosh etaD (|20o4 studied the "phase travel-time" (TT: Finsterle et alJEooll) of magneto 



acoustic waves observed to propagate between the 1700A and 1600A TRACE ultraviolet contin- 
uum pass bands and their correlation to the extrapolated chromospheric topography in an equa- 
torial coronal hole. Changes in travel-time measured inside a coronal hole correlate strongly to 
the magnetic topography present and suggest a significant change in atmospheric conditions at the 
base of the chromosphere, relative to the quiet Sun. This single observation of possible structural 
changes in the chromosphere under a coronal hole sparked a broader question: Is there a significant 
connection between the structure of the chromosphere and that observed in the solar wind? 
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Mcl ntosh & Leamon (2005), using very simple assumptions, showed that there was a signifi- 
cant connection between the inferred structure of the chromosphere and direct ACE measurements 
of solar wind conditions for a variety of chromospheric magnetic topologies. Specifically, we 
demonstrate that diagnostics of atmospheric "depth" in the chromosphere, such as the travel time 
of dispersive magneto-acoustic waves between the 1600A and 1700A TRACE UV passbands, cor- 
relate very strongly with in situ solar wind velocity (Vsw) and inversely with the ratio of ionic 
oxygen (0 7+ /0 6+ ) densities, and that a single power law relationship existed between TT and 
Vsw or 7+ /0 6+ held for fast outflows from coronal holes, slow outflows from the streamer 
belt and intermediate outflows from quiescent sun regions. The best-fit power law relationships, 
X = A(TT) B + C, are, for the wind velocity V sw , 

A = 0.053 ± 0.007; B = 4.56 ± 0.33; C = 333 ± 12 (1) 

and for the composition ratio 7+ /0 6+ , 

A = 29.9(±2.70) x 10 3 ; 

B = -7.21 ±0.23; C = 0.011 ± 0.003. ^ 

Even before it enters the corona, the solar gas is fractionated; transport through the chro- 
mosphere into the corona causes an overabunda nce of elements w ith low first ionisation potential 



such as magnesium (the so-called FIP effect). iGeiss et al.1 (|1995|) studied the variations in Mg/O 
and 7+ /0 6+ during Ulysses' 1992-93 passages through the southern High Speed Stream (i.e., 
Coronal Hole), and found that they were highly correlated. Since the Mg/O ratio is controlled by 
the FIP effect and the 7+ /0 6+ ratio reflects the coronal temperature, the correlation between the 
two points to a connection between coronal conditi ons and the structu re of the chromosphere. We 



find that Mg/O is also highly correlated with TT. IGeiss et al.l (|1995|) further showed that, as the 
Sun (and the coronal hole) rotated under Ulysses, the transition between the coronal hole and quiet 
sun was almost as steep as a step function. Thus the chromosphere and corona have a common, 
relatively sharp, boundary, separating the low-FIP from the high-FIP region in the chromosphere 
and the low-T from the high-T region in the corona. The existence of such a common boundary 
points to a causal relationship of the kind for which conditions in the chromosphere determine the 
supply of energy into the corona. 

That coronal holes give rise to fast streams of wind depleted in hot oxygen is not a new result; 
the key is that there is one continuous function describing a range of solar wind speeds derived from 
the structure of the chromosphere, just as there is a continuous energy input spectrum responsible 
for the plasma topology. 



2.1. How the footprint comes to be: 
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Energy input into the solar atmosphere 



Mci ntosh et al.l (2006b) offer an explanation of atmospheric heating and initial solar wind 
acceleration (inside and outside of coronal holes respectively) that is observationally consistent 
with a single physical mechanism being responsible for both processes. The correlations between 
different spectroscopic features observed are consistent with the energy input to the plasma via 
the action of relentless magnetoconvection-driven magnetic reconnection; the energy input to the 
plasma is modulated by the unsigned field strength on the supergranular length scale while the 
closure of the global magnetic topology determines if the injected mass and energy heats closed 
loops or is accelerated along open field lines. 

In a coronal hole, the small (spatial) scale flux elements anchored in the super-granular bound- 
ary are effectively open to interplanetary space because they have the same magnetic polarity as 
the bulk of the coronal hole. As a small, recently emerged magnetic dipole is advected to the 
boundary, the leading polarity of the flux begins to reconnect with the anchored element, creating 
a new magnetic topology in the super-granule interior. A portion of the energy released by the 
reconnection quickly begins to evaporate cool chromospheric material into the topology whereas 
the bulk of the remaining energy is released in the form of kinetic energy. 

Conversely, in the quiet Sun the magnetic flux anchored in the super-granule boundary closes 
somewhere in the vicinity of the cell with a similarly anchored, opposite polarity piece of magnetic 
flux on a nearby super-granule, such that there is very little probability that any arcade created in 
the super-granule interior can open into interplanetary space. Thus, as the reconnection progresses, 
the bulk of the released energy must result in the evaporative mass-loading of the created arcade 
and thermal heating of the plasma contained therein, since little of it can be rapidly converted into 
plasma outflow. 



Section 5.1 of iMcIntosh et al.l (I2006bl) discusses the likelihood that the energy input mecha- 
nism and its connection to the global magnetic topology are tied to Geiss' observation of the FIP 
effect. We prescribe that the fast solar wind generally originates f rom open magnetic regions (with 



sizable supergranular magnetic imbalance. IMcIntosh et al]|2006a|) and that the energy is delivered 



in a mostly non-thermal (kinetic) form with mass that originates in the well mixed chromosphere. 
Conversely, the slow solar wind originates from a largely closed magnetic topology in the quiet 
Sun, which will ensure a complex, mostly thermal, energy delivery to the plasma. It is possible 
that enhanced fractionation can take place in the plasma as it is transported through the magnetic 
topology and eventually evaporated out of the corona and into the he liosphere in a fashion that de 



pends on the exact magnetic topology of the observed plasma (e.g. JSchwadron & McComasll2003L 
and references therein). 

Many other markers of the coupling between the initial release and distribution of energy and 
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the (supergranular) magnetic flux balance of the solar plasma may have already been observed. 
For example, systematic He I 10830A line asymmetries (observed in the chromosphere at the 
solar South pole near s olar minimum) have been at tributed to a strong outflow (~8 km s _1 ) of 
chromospheric material jDupree. Penn & Jon es 1996) and occur at a time when there is a relatively 
strong unbalanced field at the pole. The magnetic flux that diffuses to the polar regions over the 
course of the solar cycle creates a long-lived "polar crown" coronal hole at solar minimum. This 
excessive imbalance in the polar magnetic field, we suspect, acts as a reservoir of kinetic energy 
for the ver y high speed winds o bserved at high heliospheric latitudes at the same phase of the 
solar cycle (IMcComas et al.lll998|) while the evolving mixture of open and closed regions and the 



resulting energy partitioning of the p lasma explain the complex heliospheric structure observed at 



other phases of the solar cycle (e.g. jMcComas et al]|2000l : ISmith et al.ll2003l) . 



The above discussion offers support to, and an energy input mechanism for, two recent solar 
wind models that relate solar wind outflow speed at 1 AU (inversely) to coronal temperatures, and 
that energy input is shared between therm al energy and wind energy in one consistent manner 
(IFiskl 120031 : ISchwadron & McComasll2003|) . Further, it justifies our hypothesis that we can use 
observations of the chromosphere as proxies for plasma conditions in the corona and solar wind, 
even if (as) the perceived direct driver of the solar wind is expansion of the million-degree corona 
high above the chromosphere. The eventual energy balance must be reflected in the travel-time 
diagnostic TRACE observations of the solar chromosphere in corona holes and quiet Sun through 
the visible changes in the plasma topography. 



3. Forecasting the Solar Wind from the Chromosphere: Method 



If we can correlate the wind speed and composition to the conditions in the solar chromo- 
sphere, can we use the correlations to predict conditions in the solar wind? In principle, the answer 
is yes, but the nature of the TRACE UV observa tions does not lend it s elf to a predictive tool: 
each of the 13 solar observations (data points) in iMcIntosh & Leamonl (|2005|) is the result of a 
1-3 hour time TRACE sequence of interleaved 1600A and 1700A images, which require careful 
co-alignment to sub-pixe l accuracy and removal of the effects of sol ar rotation. Fortunately, there 
is a striking relationship (IFinsterle et al]|2004l : IMcIntosh et al]|2004l) between the observed travel 
times and the plasma j3 (ratio of gas pressure to magnetic pressure) in the TRACE field of view at 
the (presume d) mean formation height of the TRACE 1600A passband (some 450 km above the 
photosphere; Fossum & Carlssonll2005|) . Extrapolating a co- al igned line-of-sight magnetogram 
from the Michelson Doppler Imager (MDI; IScherrer et al.lll995|) on SOHO to 450 km, computing 
the magnetic pressure (= |£?| 2 /87r) and by imposing a simple model gas pressure (|Vernazza et al 
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1981 . model VAL3C), we can compute (3 in the chromospheric region of the TRACE passbandsQ 



Fig. [His a two-dimensional histogram that shows the relationship between the observed os- 
cillation travel time and t he plasma (3 at 450 km abov e the photosphere for the suite of 13 TRACE 
observations presented in lMcIntosh & Leamonl (120051) . Using the (3-TT relationship, we can gen- 
erate a synthetic full-disk TT map from the extrapolated MDI magnetogram (3 and then use the 
power law of Eqn. (0Q) to generate a map of wind speed or other correlated in situ quantities. 
Fig.|2]shows examples the derived wind speed map, alon g with a context EUV c oronal image from 
the SOHO Extreme-Ultraviolet Imaging Telescope (EIT: iDelaboudinierdll995|) . The latter clearly 
shows the active regions and coronal holes present on the disk. 

To demonstrate our method, we compare our predictions to observations over the course of 
approximately one whole solar rotation, from July 11 to August 10, 2003. This interval was 
chosen to overlap with th e real TRACE UV observations made on July 14, 20, and 27 (see 
Mcintosh & Leamon 2005). 



Although we compute the in situ predictions for the whole Sun, we average them over an 
Earth-directed "source window," created by tracing 33 field lines back from the 2.5R & source 
surface of a Potential Field-Source Surface (PFSS) model to the photosphere. Fig. |3] shows PFSS 
field-line tracings for 4 times (separated by approximately one week) spread through our sample 
interval. The angular distance between the point the sub-terrestrial field line maps back to and the 
sub-terrestrial point ranges between zero (i.e., maps straight back) and as much as 35° away from 
disk center; the average separation is 14°. The other 32 field lines are arranged on the perimeter of 
an ellipse with semi-major axes 5° x 2.5°. These 32 field lines usually map back to a small closed 
region on the photosphere. The 4th panel of Fig. |3] shows a time when this is not the case; the 
ellipse maps to a crescent- shaped swath covering most of the northern solar hemisphere. 

To calculate th e time it takes a plasma parcel to propagate out to 1 AU, we use the method 
of ICranmerl (2004|) . invoking Lambert's W function, to solve the original solar wind model of 
Parkerl ( 1958 ). With the predicted wind speed at 1 AU as an initial condition, we integrate the 
time- of- flight backwards to the sun. Adding the time-of-flight to the MDI observation dates gives 
a time-series in situ prediction, which we compare to ACE data. Interestingly, the time difference 
between ballistic and accelerating flows arriving at 1 AU scales as V S ^J 2 . 



'The Ni I line at 6768 A used by MDI is itself formed some 200 km above the photosphere, so we are, in fact, only 
extrapolating the magnetic field 250 km. 
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4. Results and Discussion 

The lowest two panels of Fig. |4] show the results of our propagation model. For comparison 
purposes, the TRACE observations mentioned above are marked on Fig. [4] in orange. There is a 
good correlation, albeit with the predictions lagging the observations. For the whole interval, the 
peak cross-correlation coefficient of 0.47 occurs at a lag of +19. 6h for 7+ /0 6+ , and 0.32 at +24h 
for Vsw The Spearman rank-order correlation coefficients are somewhat higher: 0.63 and 0.46, 
respectively. We attribute the long lags to the unusually slow wind speeds resulting from periods 
when the wind emanates from active regions. Not unrelated is the variation in error bar size — 
smallest when the whole ellipse maps back to one point in an active region (i.e., slowest wind), 
and large (e.g., panel 4 of Fig. [3]). This is not surprising, considering that the TRACE data used to 
derive Eq. (OQ) are averaged over a 200" x 200" region and contain a range of different TT values, 
whereas the smallest PFSS footp rints are approximately 80" s quare. It is clear that a sub-field less 



than a quarter the size of that of [Mcintosh & Leamonl (|20 05 ) will derive more extreme values of 



TT and adjust the power-law relationships accordingly. 

The top panel of Fig. |4] shows the predictions for the same time interval of the LMSAL 
DeRosa-Schrijver model and the second panel those of the NOAA SEC Wang-Sheeley-Arge model. 
The LMSAL forecast is in fact 25 different predictions (red traces), from 25 slightly different ini- 
tial conditions (varying base wind speed, stream interaction strength, etc.). Similarly, the WSA 
model is uses data from three different sources as inputs. Much has been made of the intercalibra- 
tion between the Mount Wilson (MWO), Wilcox (WSO) and Kitt Peak (NSO) data used, but it is 
not germane to discuss them here. It is reassuring to see that our model, based on chromospheric 
TT, gives essentially similar results to the other predictions, where the wind speed is related to the 
super-radial expansion factor of at each point on the Source Surface. A visual inspection suggests 
that our model most closely tracks the WSO-derived predictions. For such a simple model, the 
correspondence with both observations and more complex models is both striking and encourag- 
ing. We cannot quote correlation coefficie nts for the other models for the interval studied, but the 



WSA method paper (|Arge & Pizzol 12000) cites that for the three-year period centered about the 
May 1996 solar minimum, the correlation coefficient for Vsw is 0.4 with an average fractional 
deviation of 0.15. We must be happy with 0.47 and 0.32, respectively, for one month's worth of 
predictions given the greater variability of the solar wind in the declining phase of a solar cycle 



( McComas et al.l l2000) . 



We have previously reported on the correlations that exist between in situ solar wind quan- 
tities measured by ACE and diagnostics of the chromospheric plasma topography from TRACE 
observations. The two main results of the present work are: (1) the chromospheric plasma (5 
map derived from SOHO MDI data can be used as a useful full-disk proxy for TRACE travel- 
time measurements of chromospheric structure; and (2) the plasma j3 maps can be used generate 
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accurate forward-modelling predictions of the in situ solar wind state, with correct magnitudes 
and prediction of gross features and these predictions can be updated at the cadence of the full- 
disk magnetograms. We emphasize again that this is a very simple, empirical model connecting 
observed changes in chromospheric structure to their impact on the structure of the solar wind. 
Nevertheless, as a first generation model (over a limited time interval), we obtain a reasonable 
amount of detail on the structure of the nascent solar wind. One goal of future work is clearly the 
development of a real-time model. On the other hand, perhaps we are not as truly empirical as we 
suggest: We deviate from existing wind models by predicting composition. While 7+ /0 6+ might 
not be a priority for space weather, that we can predict it with the same accuracy as Vsw affirms 
the core concepts and physics behind the chromosphere-corona coupling outlined in section [27T1 

Finally, several limitations to our model, including the using the quiet-Sun VAL3C atmo- 
spheric model everywhere and using a potential extrapolation to generate maps of (3 and TT, can 
be simply addressed by using full disk time-series observations of the chromosphere to compute 
the travel-time maps without having to perform the present limited field-of-view to full disk trans- 
form. Full disk observations of this type will be a standard "data product" of the Solar Dynamics 
Observatory. 

SWM's contribution is based upon work supported by NASA under Grants NNG05GM75G, 
issued under the Sun-Earth Connection Guest Investigator Program and NNG06GC89G, issued 
under the Solar & Heliospheric Physics Program and by the NSF through its Solar Terrestrial 
Physics Program (ATM-0541567). 
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Fig. 1 . — 2-D histogram showing the relationship betwee n plasma at 450 km and T RACE UV 
travel time, for all pixels in all of the 13 events studied by iMcIntosh & Leamonl (|2005l) . The bins 
are 0.05 wide in both TT and log 10 j3, and the color scale is in pixels per bin. The best fit curve is 
TT = -23.3 + 7.0 log 10 (/3 + 5 • 10 3 ). 



-11- 




Fig. 2. — (left) Predicted full-disk solar wind speed map, as inferred from the 2003 August 6, 
03:15 UT SOHO MDI magnetogram using the (3-TT relationship shown in Fig. [Hand the TT- 
V sw relationship of Eq. ©. (right) SOHO EIT context image from 2003 August 6, 03: 12 UT. The 
yell ow and orange contou rs show the 100 and 200 DN intensity levels in the image, respectively 
(cf. iMcIntosh et al.1 120041) . In both panels, the white dashed circle has radius 600", where the 
accuracy of the line-of-sight approximation for magnetograms degrades significantly. There are 
visible artifacts beyond this perimeter, especially towards the northeast and southwest. 
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Fig. 3. — Tracing the sub-terrestrial field lines back to the photosphere. For each MDI magne- 
togram, 33 field lines are traced back to the sun; one at the sub-solar point, and 32 around an 
ellipse 5° x 2.5°. Sometimes, the whole ellipse maps to a small region of the photosphere, other 
times it maps to regions all over the photosphere. Sometimes the sub-terrestrial point on the source 
surface is connected to the sub-terrestrial point on the photosphere, other times the footpoint is 
as much as 35° away from disk center. [See the electronic edition of the Journal for an mpeg 
animation of this figure.] 
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Fig. 4. — (lower two panels) Predicted solar wind speed and oxygen composition O + /0 + (bot- 
tom) extrapolated from inferred TT and propagated out to 1 AU. Quantities as observed by ACE 
are in green, and the three orange points are real TRACE TT observations, (upper two panels) 
Predicted wind speed from DeRosa & Schrijver (above; |http : //www, lmsal/f orecast7] ) 
and Wang, Sheeley & Arge (below; [http : / /www, sec . noaa . gov/ws/) . 



